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Abstract 

Singular value analysis can give conservative sta- 
bility margin results. Applying structure to the un- 
certainty can reduce this conservatism. This paper 
presents flight-determined stability margins for the 
X-29A lateral-directional, multiloop control system. 
These margins are compared with the predicted un- 
sealed singular values and scaled structured singular 
values. The algorithm was further evaluated with flight 
data by changing the roll-rate- to- aileron-command- 
feedback gain by ±20 percent. Minimum eigenvalues of 
the return difference matrix which bound the singular 
values are also presented. Extracting multiloop singu- 
lar values from flight data and analyzing the feedback 
gain variations validates this technique as a measure 
of robustness. This analysis can be used for near-real- 
time flight monitoring and safety testing. 
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Introduction 

Multivariable control systems have been used for 
decades; however, the methodology to evaluate the sta- 
bility margins was not developed until the last decade. 1 
Classical frequency u> analysis methods, such as Bode 
or Nyquist techniques, work well for single-input- 
single-output (SISO) systems but are inadequate for 
multiple-input-multiple-output (MIMO) control sys- 
tems. Classical methods do not allow for simultane- 
ous variations of phase and gain in all of the feedback 
paths. 1-3 Recently, singular value a norms of the re- 
turn difference matrix (RDM) have been considered a 
measure of the system stability margin for multiloop 
feedback control systems. 1,2,4 However, singular value 
norms of systems with unstructured uncertainty can 
be overly conservative, and a control system designer 
could interpret the results as unsatisfactory when, in 
fact, the system is robust. 3 A method for relieving 
the excessive conservatism is derived by structuring the 
uncertainties. 2,3,5 A control system design is robust 
when it can perform well for substantial variations in 
plant dynamics. 

Although mathematically sound, singular value 
methods require substantial evaluation using multiloop 
control systems before these methods can be accepted 
as adequate for determining the robustness of control 
system design. Some insight into application of experi- 
mentally determined singular values of a multiloop flut- 
ter suppression control system for a wind-tunnel aeroe- 
lastic model is presented in Ref. 6. In Ref. 6, singular 
values of the RDM at the input and output locations 
were used successfully to evaluate the performance of 
the control system. Reference 6 also gives computa- 
tions for the open-loop transfer matrix GH from open- 
loop control system wind-tunnel operations. This ref- 
erence concentrated on wind-tunnel aeroelastic control 
performance. Evaluation is also needed for unsealed 


2 


singular value (USV) and scaled structured singular 
value (SSV) methods for analyzing multiloop control 
systems of manned aircraft in flight. In addition, theo- 
retical and flight-deter mined stability margins need to 
be compared. 

To evaluate the stability of an in-flight multiloop sys- 
tem, the NASA Dry den Flight Research Facility con- 
ducted a flight test program on a MIMO flight control 
system using the X-29A aircraft (Grumman Aerospace 
Corporation, Bethpage, New York). The X-29A is an 
experimental aircraft designed to test the integration 
of several technologies. This aircraft was ideal for this 
study because the lateral-directional control system is 
multiloop with the capability of controlled excitation 
and measured response for MIMO analysis. 7 

This paper presents an extension of the evaluation 
method reported in Ref. 6 by determining singular val- 
ues from flight data and comparing the values with the 
predicted or analytical USV and SSV, (Jjjsv an< ^ 2-SSV* 
Flight results of the minimum eigenvalues A obtained 
from the RDM and comparisons with the analytical 
minimum eigenvalues A analy are P rcsen ted. To 
determine the sensitivity of the MIMO stability margin 
algorithm, a lateral-directional X-29A control system 
feedback gain was changed during flight. Several flight 
conditions were flown and analyzed, but for the sake of 
brevity, one representative case is reported here. The 
flight test results are for three control law gain settings 
(80- , 100-, and 120-percent K 2) at a flight condition 
of Mach number M 0.7 and an altitude h of 30,000 ft. 
The technique described in this paper can be used to 
obtain near-real-time stability margins during flight for 
multiloop control systems. 

Test Vehicle Description 

The X-29A is an experimental aircraft designed to 
test the feasibility of integrating several modern tech- 
nologies into a highly maneuverable aircraft. This 
single-seat aircraft is powered by an F404-GE-400 en- 
gine (General Electric, Lynn, Massachusetts) with a 
thin supercritical 30° forward-swept wing and close- 
coupled canards. The wing contains full-span, double- 
hinged, trailing-edge flaperons which also provide vari- 
able camber. The wing structure includes aeroelasti- 
cally tailored graphite-epoxy covers designed to provide 
stiffness to overcome the structural divergence prob- 
lems associated with forward-swept wings. Figure 1 
shows the X-29A aircraft. For a more complete de- 
scription of the aircraft, see Refs. 7, 8, and 9. 

Longitudinal control of the aircraft is achieved with 
canard, symmetric flaperons, and strake surfaces. The 
aircraft was designed with a high degree of longitu- 
dinal static instability, 35 percent at subsonic speeds. 
This instability requires the vehicle to be continuously 


stabilized by a fly-by-wire control system. Lateral 
directional motion is controlled by a conventional 
rudder and differential flaperon deflection. The X-29A 
lateral directional control system is multiloop and in- 
cludes aileron-to-rudder and rudder-to-aileron inter- 
connects for turn coordination. 

Flight Control System Description 

The X-29A aircraft has a triplex digital flight con- 
trol system with an analog backup for each channel. 
The digital control law outputs are computed at 40 Hz 
in dual- processor flight computers. The primary task 
of the longitudinal control system is to stabilize the 
motion of the aircraft. Since the longitudinal flight 
control system is an SISO system, it is not evalu- 
ated in this report. The X-29A high angle- of- attack 
a, lateral-directional control system is multiloop. Fig- 
ure 2 shows the major digital and dynamic elements. 
There are two pilot inputs: lateral stick and rudder 
pedal. The four aircraft feedback signals are roll rate 
p, yaw rate r, roll attitude </>, and lateral accelera- 
tion Ny. As Fig. 2 shows, seven paths are summed to 
generate the aileron command, including the ground- 
generated signal (GGS) command. Also, seven paths 
are summed for the rudder command, including the 
GGS command. The multiloop lateral directional con- 
trol and GGS systems, which can be summed to the 
actuated commands, make the X-29A ideal for MIMO 
stability analysis. 

Roll-rate- to- aileron- feedback gain K 2 is a function 
of flight condition. The A 2 was the only gain that 
could be changed in flight that had any impact on the 
closed-loop stability margins. The K 2 shown in Fig. 2 
could only be changed by ±20 percent of the nominal 
value to determine the sensitivity of the singular value 
algorithm. The pilot could change the control system 
software K2 in flight by selecting a set of predetermined 
gains. The X-29A dynamics could also be excited by 
remote signals from a ground-based computer using the 
GGS system. For this study, the important control 
system is located at the input node where the signals 
are summed before being sent to the actuators (f ig. 2). 
The aileron-command-to-actuator 6 a cmd and rudder- 
command-to-actuator S r cmd signals are the resultant 
summed commands from the pilot, GGS, and control 
system feedback paths. 

Ground-Generated Signal System 
Description 

The X-29A flight control system was designed with 
the ability to remotely excite the aircraft dynamics 
with sweeps, step inputs, and doublets using a ground- 
based computer linked to the control system by teleme- 
try transmission. The GGS system is used to repeat 
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maneuvers, provide for independent control surface 
inputs, and produce well-defined inputs for frequency 
response calculations. 9 Figure 2 shows that the GGS 
commands can be included in the aileron or rudder 
paths. For this analysis, the GGS consisted of a sinu- 
soidal frequency sweep from 40.0 to 0.1 rad/sec. Refer 
t° Ref. 9 for a more detailed discussion of the GGS 
system. 

Data Acquisition System Description 

The GGS, aileron, and rudder commands to the ac- 
tuator are the only digital flight control system signals 
needed for this MI MO analysis. These signals were 
integrated with other data from a data bus and down- 
linked to the ground at 40 samples/sec. Figure 2 shows 
the GGS, aileron, and rudder commands to the ac- 
tuator. The measured data were relatively noise free 
because the pilot inputs and feedback paths included 
antialiasing filters. For more information about the in- 
strumentation, see Ref. 10. 

Analytical Methods 

Two methods of determining the analytical multi- 
variable stability margins are discussed in this section. 
The two analytically determined robustness methods 
are referred to as USV and SSV. 3,5 The USV method 
can be too conservative, and the results could be in- 
terpreted as unsatisfactory stability margins when, in 
fact, these margins are adequate. The SSV method 
decreases the conservative nature of singular values. 
Structuring or scaling was not added to the flight- 
determined singular values because the flight data con- 
tains the complete system, including nonlinear and 
high-order system dynamic effects. The following sub- 
sections provide background information on stability 
margins using eigenvalues, flight data, and singular val- 
ues as well as analytical structured and unsealed sin- 
gular values. 

Stability Margins Using Singular Values 

Singular value analysis has been the focus of con- 
siderable interest in the controls discipline. 1 ~ 6 > n This 
analysis provides a way to determine how much uncer- 
tainty can be tolerated before a multivariable system 
becomes unstable and assumes that the system is ini- 
tially stable. The singular values of the RDM arc used 
to measure the stability margins. 1 ' 4 The RDM at the 
input node is [I + HG]. At the output node, the RDM 
is [I + GH]. As the minimum singular value a of the 
input or output RDM approaches zero, the system be- 
comes increasingly less stable. 

lo analyze the RDM, a review of the general singular 
value analysis is helpful. Let A be a general [n xn] ma- 
trix, then the minimum and maximum singular values 
a are 


a[Ax] = min |j Ax |j = v%nin[A*A] 

Mll = i (i) 

[Ax] = max || Ax || = v / A max [A*A] 

II x || - 1 (2) 

where |j Ax || is the Euclidean norm; A [A* A] are the 
eigenvalues; and A* is the conjugate transpose of ma- 
trix A. Note that the vector Ax depends on the units 
of the control system output vector x variables; there- 
fore, the singular values also depend on the scale units 
of x variables. 2 Other useful properties of singular val- 
ues are as follows: 2,3 

dA-^lMA] (3) 

a[A] < | A[A] | < a[ A] (4) 

Equation (4) states that the magnitude of the eigenval- 
ues of A is bounded by the minimum and maximum 
singular values of A. 

Figure 3 shows a typical control system. The G(s) 
is the plant transfer function matrix, and H(s) is the 
control system transfer function matrix. Bold-faced 

variables, such as G and H, are multiple path input 

and output matrices. The Lj and Lo represent the 
input and output perturbation matrices. The MIMO 
systems have cross-feed interactions which cause the 
locations of perturbations to impact the singular value 
results. In classical SISO systems, locations of the dis- 
turbances are unimportant. Ignoring the perturbation 
matrices for now, the relationships of the RDM for the 
input node [I + HG] and the output node [I + GH] 
are usually different. Therefore, the singular values are 
different because matrix multiplication is not commu- 
tative, i.e., GH ^ HG. 

If the multiloop system is stable when unperturbed, 
then a sufficient but not necessary condition for the sys- 
tem to remain stable when Lj is perturbed is that 2,11,12 

d(Li -1 - I) < a [ I + H(;w)G(jw)] (5) 

When analyzing the input node, the output perturba- 
tion matrix is set to identity, Lo = I. Likewise, the in- 
put perturbation matrix Lj is set to identity when ana- 
lyzing the output node. The right-hand side of Eq. (5), 
a[\ + HG], is the minimum singular value of the input 
node RDM and is used as the basis of this robustness 
analysis. Analysis of singular values will work for SISO 
systems for determining stability margins. 

Analytical Unsealed Singular Values 

A fully populated perturbation matrix Lj produces 
the most conservative stability margins. The robust- 
ness of a system with a fully populated Lj can give an 
unrealistic measure of the stability margins. 1-3 Often, 
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these perturbations can be adequately defined using 
the diagonal elements of L*. When the Lj matrix is 
diagonal, it takes on the structure 

1^ = Diag[fc,e^ 1 , k 2 c . . . k n c**') (6) 

where k x is the uncertainty gain element, and <f> x is the 
uncertainty phase element of the L matrix. Also, k x 
and <f> x (z = l.n) may change independently and simul- 
taneously with respect to one another in each control 
loop. This method allows for the robustness analysis 
of multivariable control systems. 

Now that the ground work has been laid, the sys- 
tem stability boundary can be defined by testing the 
criterion 

£[I + HGLi(jw)] = 0 (7) 

as a function of frequency. This expression is the 
neutral stability boundary that produces a pure res- 
onance at lj . This relationship means that a pole of 
the closed-loop system on the imaginary axis exists. In 
addition, the minimum singular value of the RDM is 
zero. Whether the closed- loop system is stable or un- 
stable, the singular values are always nonnegative as 
Eq. (1) indicates. The singular values are magnitudes 
of a transformation that are always positive and define 
the distance to the neutral stability boundary. 3 The 
reference for neutral stability is a minimum singular 
value of 0.0, and the stability margins are related to the 
magnitude of minimum singular value. 1 he greater the 
magnitude of the minimum singular value, the greater 
the robustness. 

One goal of this research was to extend the well- 
understood, classical stability margin methods from 
SISO to Ml MO controllers. Equation (7) represents 
the information needed to determine the stability mar- 
gins of a system. However in this form, this equation 
does not relate the magnitude of the minimum singular 
value to traditional phase and gain margins. The rela- 
tionship that expresses the gain and phase as a direct 
function of either the input or output perturbations 
is shown in Eq. (8) and derived in Ref. 4, where L is 
either L| or L 0 

cr(L -1 — I) = y ^1 — — ^ — cos(<jS>n ) ) W 

Figure 4 shows a graphical representation of Eq. (8) 
and presents minimum singular value as a function 
of phase and gain margins. To determine the stabil- 
ity margins or nearness to instability, compute a[I T 
HG(jtj)] for various frequencies. The minimum singu- 
lar value plot traces the nearness to singularity of the 
RDM and, thus, the system robustness to perturba- 
tions as a function of frequency. It is not necessary 
to know the perturbation matrix L to compute the 


minimum singular value plot. For more information 
regarding the stability boundary and its relationship 
to the perturbation matrix L, see Ref. 3. 

Analytical Sealed Structured Singular Values 

Because the singular values of a fully populated L- 
matrix arc conservative, it may be misleading to apply 
these values to a control system. Such values may indi- 
cate that the system is not robust. The Analytical Un- 
sealed Singular Values subsection structured the uncer- 
tainty bv making the L- matrix diagonal (see Eq. (6))- 
The USV method reduces the conservatism; however, 
scaling the system further reduces the conservatism. 
Scaling the system provides a way to reduce the con- 
servatism and still maintain realistic margins. 1,3 ’ 0 I his 
reduction is accomplished by including a diagonal scal- 
ing matrix D in the RDM expression. The scaling 
should be chosen so jus to maximize the minimum singu- 
lar value across the frequency range. Therefore, D(u/) 
is a function of frequency, and the algorithms are easily 
implemented. 13 

a [I + D(jw){H (^)G0^)}D(^) -1 ] < SSV (9) 

Equation (9) represents the singular values as a func- 
tion of frequency in the presence of the matrix D(cj), 
which reduces the sensitivity of cross-feed perturba- 
tions. Implementation of the D matrix minimizes the 
conservative nature of singular values of multiloop ro- 
bustness predictions. Equation (9) eliminates the con- 
servatism for control systems where the dimension of 
HG is three or less. 5,11 Scaled singular values arc the 
same as SSV and are used only for the analytical por- 
tion of the analysis. The predicted singular values of 
the RDM, [I + HG] as well as [I + DHGD" 1 ], can 
be determined from linear lrequency analysis. 

Eigenvalue Analysis 

As shown in the Stability Margins Using Singular 
Values subsection, the node point location of the anal- 
ysis can influence system stability. In addition, the 
singular values of the RDM can give a conservative 
measure of the closed-loop stability margins. Note that 
system singular values will always be upper bounded 
by the eigenvalues of the RDM, see Eq. (4). 1 he eigen- 
values of the RDM are identical at any location in the 
control system 6 and can be expressed as 

| A(I + HG) | = | A(I + GH) | (10) 

input node output node 

Equation (10) represents the upper limit of the mini- 
mum singular value and, therefore, the upper bounds 
of the stability margins of a multiloop control sys- 
tem. The eigenvalues of the RDM are not the same as 
those of the closed-loop system. As a reminder, closed- 
loop eigenvalues must be stable before applying this 
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stability margin analysis technique because minimum 
singular values only indicate how far from neutral sta- 
bility the system is either stable or unstable. 

Singular Value Analysis of Flight Data 

The flight singular values and eigenvalues need to 
be determined by using frequency response techniques. 
This section describes the methodology used to deter- 
mine singular values from flight test data. Transfer 
function equation development, eigenvalues, and sin- 
gular values of the RDM are included. 

A complex frequency response of a system can be es- 
timated from the autospectrum and cross spectrum of 
the input and output time history variables by trans- 
forming these time domain responses to the frequency 
domain using fast Fourier transforms (FFT’s). The 
controller input-to-output, (u(*)}-to-{x(*)}, transfer 
matrix (fig. 3) Xu is defined as follows: 

N 

{Xut?'w)},j =53( Sx > u i(J w ))i{Su i u i e-w),}- 1 (11) 
1-1 

where, S X u is the cross spectrum of the input u and 
output x. The Suu is the autospectrum of the input, 
and N is the number of time history arrays. The data 
were loaded into arrays of 2048 points. A raised cosine 
smoothing window was used to process the time history 
data obtained from flight. See Ref. 8 for more details of 
the FFT procedures. The transfer functions produced 
by GGS excitation yielded good results. 

Response matrix Xu can be used to construct the 
RDM [I + HG], The following development shows the 
relationship of X u to the RDM. Matrix terms HG or 
are called the loop gain matrices. Closed-loop sys- 
tem robustness is determined from the minimum singu- 
lar values of the RDM at cither the plant input node, 
2-1 1 + HG], or the output node, a [I + GH], As the 
minimum singular values approach zero, the stability 
margin goes to zero. Perturbation matrices are not 
required for flight-determined singular value analysis 
because variations are inherent in the system dynam- 
ics. The following closed-loop relationships can be de- 
veloped from Fig. 3 if the perturbation matrices are 
ignored: 

e = u-x (12) 

x = HGe = HG[u - x] (13) 

The complex frequency response of the open-loop 
transfer function can be estimated from the autospec- 
trum and cross spectrum of u and x. 

Sxu = HG[S UU - S xu ] (14) 

This relationship can be postmultiplied by the au- 
tospectrum inverse Suu -1 - 

SxuSuu' 1 = HG[S uu Suu 1 - SxuSuu' 1 ] (15) 


Combining Eqs. (11) and (15) produces 

X u = HG[I - X u ] (16) 

Therefore, the loop gain matrix as a function of fre- 
quency is 

HG(jcu) = X U (^)[I - XuO'w)] -1 (17) 

The response matrix X u (jw) is square and has the 
dimensions equal to the number of inputs. The ex- 
pressions have now been developed to determine the 
flight test singular values for the stability margins of 
multiloop control systems. The spectral relationships 
of this section can be rapidly evaluated using FFT’s 
which make it possible to determine, during flight, the 
near-real-time stability margins of multiloop control 
systems. 

Flight Test Procedure 

The maneuvers flown for the multivariable stabil- 
ity margin analysis were designed to excite the mo- 
tion of the lateral-directional axis. The pilot would 
stabilize the aircraft at the desired flight condition. 
Then, a GGS frequency sweep (Fig. 5) would be com- 
manded to the roll axis to excite the X-29A dynam- 
ics, The GGS maneuver was complete after approx- 
imately 45 sec, and the pilot would reestablish the 
initial flight condition. Next, the same GGS frequency- 
sweep signal would be commanded to the yaw axis 
for approximately 45 sec to complete the needed in- 
puts for the transfer function estimation. The GGS 
frequency-sweep signal started at 40 rad/sec and fin- 
ished at 0.1 rad/sec. This type of signal sweep helped 
to maintain the initial flight condition without re- 
quiring pilot corrections. Ground-generated as well 
as aileron- and rudder-commands-to-actuator signals 
were recorded for the MIMO analysis at 40 samples/sec 
during the maneuvers. For the three K2 values, the pi- 
lot would dial in the appropriate settings on the control 
system panel, and the GGS maneuvers were repeated. 
The X-29A flight control system with its interconnects 
will not allow the pilot to generate independent con- 
troller excitation. Such excitation is essential for multi- 
loop frequency analysis. As a result, pilot-commanded 
maneuvers were not used for this analysis. 

Results and Discussion 

This section presents flight results of singular values 
as a function of frequency and compares these results 
with the analytical USV and SSV. The minimum eigen- 
values of the RDM and an evaluation of how the sin- 
gular value analysis algorithm performed when a single 
gain in the feedback path was changed by ±20 percent 
are also presented. As mentioned in the Flight Control 
System Description section, K 2 was the only gain that 


6 



could be changed easily in flight that had any impact 
on the stability margins. In addition, the flight condi- 
tions were M = 0.7 and h = 30,000 ft with K 2 of 80, 
100, and 120 percent of the nominal value. 

The ground- generated input signals to the roll and 
yaw axes resulted in the response shown in Figs. 6 and 
7. The robustness analysis for the X-29A was at the 
input node [I + HG] and had two inputs (roll and 
yaw GGS commands) and two outputs (aileron and 
rudder commands to actuator). Therefore, the RDM 
dimension will be a 2 x 2 matrix that is a function of 
frequency. 

Figure 8 shows the flight-determined input node min- 
imum singular values, q _[ I T HG], as a function of fre- 
quency with the nominal K2 as well as analytical struc- 
tured and unsealed singular values. The plot shows 
that good agreement exists between the flight and ana- 
lytical data. The analytical SSV’s tend to agree slightly 
better with the flight data than the analytical USV’s. 
This result is consistent with the theory. Figure 9 
shows the scaling or structuring used on the analytical 
system for the nominal K 2 maneuver. The dimension 
of the diagonal D is 2 x 2, which is the same size as 
HG. The scaling algorithm set the D 22 element to 1.0 
and left the D n element free to maximize Fq. (9). As 
shown in Fig. 9, Du is the same order of magnitude as 
T)' 2'2 Therefore, the scaling will not have as large an 
impact on the results as it would if the scaling were to 
differ by several orders of magnitudes. Another point 
to note is that as the D matrix approaches unity I, the 
SSV’s become the same as the USV’s. Comparing the 
minimum singular values in Fig. 8 using scaling D ele- 
ments in Fig. 9 shows that the modeling of the X-29A 
analytical system was good. 

Obtaining good low-frequency results is difficult 
when the flight maneuver is time limited (45 sec for 
this analysis) because several periods at a given fre- 
quency may be needed to accurately define the trans- 
fer function. The X-29A flight data transfer functions 
were truncated below 0.9 rad/sec because of increas- 
ingly poor transfer function definition. However, the 
important location on a singular value plot is the global 
minimum singular value, which is 0.72 at 8.0 rad/sec 
(Fig. 8), because this represents the worst closed-loop 
stability margin. The universal phase and gain margin 
plot of Fig. 4 is required to relate the g_ = 0.72 to a 
stability margin. From Fig. 4, a singular value of 0.72 
corresponds to gain margins of —4.8 and 11.5 dB and 
to a phase margin of ±41°. The dashed lines high- 
light the area of interest in this figure. These singular 
values thus imply that the gain in both paths can be 
increased by 11.5 dB or reduced by 4.8 dB simultane- 
ously before the system becomes unstable. Similarly, 
the phase in both paths can be changed by ±41° before 


the system becomes unstable. Flight and analytical 
scaled global minimum singular values are close to each 
other. The analytical unsealed global minimum singu- 
lar value is approximately 0.65 at 8 rad/sec, which cor- 
responds to a gain margin between -4 and 8.5 dB and 
a phase margin of ±35°. The SSV method matched 
the peak and valley of the flight minimum singular val- 
ues better than the USV method. As expected, the 
analytical USV stability margins were lower or more 
conservative than the analytical SSV margins. 

Figure 10 shows the minimum eigenvalue as a func- 
tion of frequency of the RDM [I + HG]. Minimum an- 
alytical eigenvalues compare very well with the flight- 
determined RDM eigenvalues. The eigenvalues are in- 
variant under scale changes (as mentioned in the Eigen- 
value Analysis subsection (Eq. (10))). The minimum 
eigenvalue is the upper bound of the minimum singu- 
lar value of the RDM, as expressed in Eq. (4). Fig- 
ure 11 combines the flight-determined minimum eigen- 
value from Fig. 10 and the minimum singular value 
from Fig. 8. The singular value curves are equal to or 
below the eigenvalue curve which agrees with Eq. (4) 
[a < A]. Analytical USV’s are not presented. 

To investigate the ability of the algorithm to detect a 
change in stability margins, K 2 was changed in flight 
by ±20 percent, and the same flight maneuvers were 
performed again. Figure 12 shows the results for K2 
of 80 percent of the nominal value flown at M = 0.7 
and h = 30,000 ft. Figure 13 shows the results for K 2 
of 120 percent of the nominal value. The analytical 
USV’s are generally lower than the flight singular val- 
ues and are always less than or equal to the SSV’s. The 
SSV curve compares well with the in-flight values at the 
lower frequencies. In addition, the USV and SSV are 
in close agreement at the minimum singular values lo- 
cation. For comparison, Fig. 14 shows the three flight- 
determined singular values for the three K 2 settings: 
80, 100, and 120 percent. The algorithm detected the 
multiloop stability margin change caused by the single 
K2 change of ±20 percent. 

Since the analytical results matched the flight results 
very well, it can be concluded that the flight minimum 
singular value can be used as a measure of multiloop 
stability margins. Although the near-real-time analy- 
sis was not done in this report, the time used to gener- 
ate flight-determined minimum singular value was less 
than 30 sec and is considered insignificant. The size 
of the RDM [I + HG] was 2 x 2 x 1024 frequen- 
cies. These computation times are small enough to sup- 
port near-real-time, multiloop stability margin analy- 
sis. The near- real- time capability would minimize the 
time required for envelope expansion of aircraft with 
multiloop control systems. 
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Concluding Remarks 

Multiloop stability margins were determined for the 
X-29A aircraft from flight data using the methods 
presented in this paper. The flight results compared 
well with predicted stability margins. Analytical sta- 
bility robustness was determined using unsealed and 
scaled structured singular value analyses. The flight- 
determined singular values were determined using the 
closed-loop frequency responses. Data analysis com- 
paring predictions of both methods showed good cor- 
relation. However, the scaled structured singular value 
method matched the flight minimum singular value 
better at the lower frequencies. The predicted un- 
sealed singular value mini mums were always conserva- 
tive compared with the scaled structured singular val- 
ues throughout the entire frequency range. Predicted 
and flight-determined minimum eigenvalues of the re- 
turn difference matrix were also presented. Sensitivity 
of the algorithm was evaluated by changing a feedback 
gain by ±20 percent, and the stability margins were 
compared with the nominal gain results. The analy- 
sis method is suitable for detecting changes in stability 
margins. 

Extracting multiloop singular values from flight data 
and comparing the information with prediction val- 
idates the application of the technique as a relative 
measure of robustness. This comparison increases the 
confidence of using singular values for stability assess- 
ments of multiloop control systems. The technique de- 
scribed in this report can be used on any multiloop 
control system. Also, this technique extends the single- 
loop gain and phase margin concepts to multiloop 
systems. The methodology could be implemented in 
near real time for flight-monitoring and safety require- 
ments. Near-real-time capability would minimize the 
time required for envelope expansion of aircraft with 
multiloop control systems. 
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Fig. 1 The X-29A aircraft. 
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Fig. 2 Lateral-directional control system. 
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Fig. 3 Control system diagram with disturbance ma- 
trices at the input and output nodes. 



Fig. 4 Universal diagram for multiloop gain-phase 
margin evaluation. 
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Fig. 5 The ground-generated signal frequency-sweep input for roll and yaw axes 
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Ftg 6 Roll axis frequency-sweep time histories caused by ground-generated signals telemetered from the ground 

to the aderon-path input for a Mach number of 0.7, an altitude of 30,000 ft, and a roll-rate-to-aileron-feedback 
gain of 100 percent. 
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Fig. 7 Yaw axis frequency-sweep time histories caused by ground- generated signals telemetered from the ground 
to the rudder-path input "for a Mach number of 0.7, an altitude of 30,000 ft, and a roll- rate- to- aileron- f eedback 
gain of 100 percent. 
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Fig. 8 Flight and predicted minimum singular values as a function of frequency for a Mach number of 0./, an 
altitude of 30.000 ft, and a roll-rate-to-aileron-feedback gain of 100 percent. 
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Diagonal scaling matrix D(w) applied to the input return difference matrix for a Mach number 
e of 30,000 ft, and a roll-ratc-to-aileron-feedback gain of 100 percent. 
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I 4 ^ Flight and predicted minimum eigenvalues of [I + HG] for a Mach number of 0.7, an altitude of 30,000 ft 
and a roll-rate- to-aileron- feedback gain of 100 percent. 
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Fig. 13 Flight and predicted minimum singular values of input node [I + HG] for a Mach number of 0.7, an 
altitude of 30,000 ft, and a roll- rate- to- aileron- feedback gain of 120 percent. 



Fig. 14 Flight minimum singular values as a function of gain settings for a Mach number of 0.7, an altitude of 
30,000 ft, and a roll-rate-to-aileron-feedback gain of 120, 100, and 80 percent. 



